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ABSTRACT 


Binary  sequences  find  increasing  use  in  electrical 
engineering  applications  of  ranging,  time  measurement  and 
communications.  A property  of  interest  in  these  applications 
is  the  autocorrelation  function  of  the  binary  sequence  or 
pair  of  sequences.  Of  the  2n  possible  sequences  of  length  n, 
only  a few  have  usable  autocorrelation  functions.  There  is, 
to  date,  no  procedure  known  which  will  provide  the  sequence 
having  a specific  autocorrelation  function,  except  for  very 
particular  cases  . 

In  this  report,  known  properties  of  complementary  sequences 
are  reviewed.  Almost  complementary  sequences  are  defined  and 
the  procedure  to  obtain  them  is  outlined.  A formula  is 
derived  for  the  number  of  different  autocorrelation  functions 
of  the  2n  possible  sequences  of  length  n bits.  A computer 
search  is  implemented  with  the  objective  of  discovering  se- 
quences with  desirable  autocorrelation  functions. 
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I.  INTRODUCTION 


This  study  is  concerned  with  binary  sequences  and  their 
autocorrelation  functions.  The  objective  is  to  obtain  auto- 
correlation functions  with  sidelobe  levels  less  than  or  equal 
to  predetermined  values.  Sequences  or  groups  of  sequences 
which  can  provide  this  property  are  very  attractive  for  use 
in  systems  whose  performance  depends  on  the  autocorrelation 
function  magnitude.  Such  systems  are  used  in  communications, 
ranging  and  time  measurement.  The  problem  is  to  find  these 
sequences  or  "good"  codes  and  the  rules  to  construct  them,  if 
such  rules  exist. 

A class  of  such  codes  are  the  complementary  sequences. 

They  are  pairs  of  sequences  with  the  characteristic  that  the 
sum  of  their  autocorrelation  functions  is  a waveform  that  has 
no  sidelobes.  These  sequences  were  first  considered  by 
M.  J.  Golay  [Ref.  l]  and  further  investigated  by  S.  Jauregui 
[Ref.  2].  They  are  used  here  to  develop  a technique  for  con- 
structing  another  class  of  "good"  codes,  the  almost  comple- 
mentary sequences. 

A.  PLAN  OF  THE  RESEARCH 

The  efforts  to  solve  the  problem  follow  two  directions: 

(1)  Experimentation  with  the  binary  sequences  and 
their  properties,  to  find  the  rules  which  give  desirable  auto- 
correlation functions. 
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(2)  Computer  search  of  binary  sequences  of  several 


lengths  to  find  the  ones  with  small  autocorrelation  sidelobes. 
As  a result  the  following  were  achieved: 

(1)  Discovery  of  the  almost  complementary  sequences. 
These  are  pairs  of  binary  sequences  which  can  be  constructed 
using  complementary  sequences.  Their  autocorrelation  functions 
when  added  have  sidelobes  of  predetermined  magnitude,  polarity 
and  position. 

(2)  Discovery  of  codes  of  lengths  £20  which  have  auto- 
correlation functions  with  sidelobes  less  or  equal  to  one. 

For  example,  for  length  n = 20 , only  three  such  codes  were 
found  and  for  n = 15,  none.  The  computer  programs  used  here 
can  be  used  to  select  codes  with  any  sidelobe  levels. 

Other  results  of  interest  are: 

(1)  Derivation  of  the  formula  for  the  number  of  the 
different  autocorrelation  functions  of  all  sequences  of 
length  n. 

(2)  Construction  of  computer  programs  which  can  be 
used  in  other  cases  as  well.  For  example,  an  algorithm  for 

, the  automatic  production  of  all  possible  binary  numbers  of 

length  n was  devised.  This  algorithm  can  be  used  to  select 
codes  having  certain  properties,  such  as  codes  with  a fixed 
number  of  ones  and  zeros . 

3.  PLAN  OF  THE  REPORT 

Chapter  II  provides  the  necessary  background  by  giving 
the  definitions  and  basic  properties  of  the  binary  sequence 
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and  autocorrelation  function.  The  formation  of  the  autocorre- 
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lation  function  Rv(t)  of  a two  level  function  v(t)  is  developed. 
Matched  filters  and  their  realizations  are  also  discussed. 

Complementary  sequences  and  their  basic  properties  are 
reviewed  in  Chapter  III,  to  form  the  basis  for  the  material  on 
almost  complementary  sequences . 

In  Chapter  IV,  the  almost  complementary  sequences  are 
defined,  and  rules  for  their  construction  are  given. 

The  formula  for  the  number  of  different  autocorrelation 
functions  in  all  sequences  of  length  n is  developed  in 
Chapter  V . 

Computer  programs  and  their  algorithms  are  discussed  next 
in  Chapter  VI . 

Chapter  VII  gives  a few  applications. 


J 
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II.  BACKGROUND 


In  this  section  some  basic  concepts  are  discussed  and 
definitions  given  as  background  material. 

A.  BINARY  SEQUENCE 

A binary  sequence  is  a list  of  elements  each  of  which  can 
have  one  of  two  distinct  values.  These  values  are  usually 
represented  either  by  +1  and  -1  or  by  1 and  0 . Often  when 
the  +1  and  -1  convention  is  used,  the  ones  are  omitted  and 
only  + and  - are  written. 

For  example,  sequence  A can  be  written: 

A = ++-+  or  A = 1101  or  A = +1+1-1+1 

The  number  of  elements  in  a sequence  is  the  length  denoted 
here  by  n.  For  the  above  example,  n = 4. 

The  voltage  equivalent  v(t)  of  a binary  sequence  is  a 
time  waveform  where  1 is  represented  by  a voltage  level  +V 
and  0 by  a voltage  level  -V. 

For  the  sequence  A = 1101,  v(t)  is  given  in  Fig.  1 where 
t is  the  bit  duration. 
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The  way  to  find  the  autocorrelation  function  Rv(t)  of  a 
digital  sequence  v(t),  is  to  "slide"  the  sequence  past  itself 
to  the  right  and  left  and  at  each  position  form  the  product 
of  the  sequence  and  its  shifted  replica.  Then  the  area  of 


the  product  waveform  is  taken  and  this  gives  the  autocorrelation 
of  the  sequence  at  this  position. 


The  procedure  is  illustrated  in  Fig.  2,  by  showing  the 
complete  steps  for  two  "shift"  positions  x = 0 and  x = e. 

It  can  be  seen  from  the  equation  for  Ry(x)  and  from  Fig.  2 
that  when  v(t)  is  a piece-wise  constant  function,  Rv(x)  will 
be  a piece-wise  linear.  The  linear  segments  terminate  at 
multiples  of  e,  a bit  duration. 

The  autocorrelation  is  an  even  function  which  has  its 

maximum  value  at  x =0.  R^tO)  gives  the  level  of  the  main 

lobe.  Secondary  maxima  are  the  sidelobe  levels.  For  the 

2 

example  of  Fig.  2,  the  mainlobe  level  is  Rv(0)  = 4V  s. 

Sidelobe  levels  are  Rv(e)  = -V^e  and  Rv(3e)  = +V^e  . 

Actually  the  shape  of  the  autocorrelation  function  of  a 
sequence  v(t)  is  obtained  easier  by  letting  v = 1 and  e = 1. 

For  example,  to  form  the  autocorrelation  function  of  the 
sequence  A = 1101,  the  sequence  is  written  and  its  delayed 
version  is  placed  beneath.  For  example,  x = 0 gives 

r 

noi ; 

1101' 

In  each  position  the  elements  of  these  two  similar  se- 
quences are  compared.  If  they  are  the  same  (both  zeros  or 
ones)  they  form  an  agreement;  if  not  (one  zero,  one  one)  they 
form  a disagreement.  The  number  of  disagreements  is  subtracted 
from  the  number  of  agreements  and  the  result  is  proportional 
to  the  autocorrelation  function  at  this  position.  Here  there 
are  only  four  agreements  and  the  result  is  4. 

L6 


j 


V ( t ) 
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v(t)v(t) 


£ 2e  3e  4e 
(a)  t = 0,  Rv(0)  = 4V2e 
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(b)  t = £,  Rv(e)  = -V2£ 


Now  a shift  is  made  as  follows 


1101; 

110 ' 1 

and  by  the  same  method  there  are  one  agreement  and  two  dis- 
agreements, so  the  result  is  -1. 

Similarly,  the  next  position  is 


1101 


which  gives  0 . 

The  next  position 


1101 


gives  1 . 


And  finally,  the  last  position 


1101 


1101 


always  gives  0 


Since  the  autocorrelation  is  an  even  function,  Rv(t)  = 
Rv(-t).  If  t is  considered  to  provide  a shift  to  the  future, 
-t  is  a shift  to  tl\ie  past. 

It  is  not  necessary  for  the  values  given  by  shifting  to 
the  left  to  be  written,  because  they  are  the  same  with  the 
ones  resulting  from  shifting  to  the  right.  So,  by  convention 
the  autocorrelation  of  the  sequence  A can  be  written 

Rv(t)  = 4, -1,0, 1,0. 

This  convention  will  be  followed  in  the  rest  of  this  report. 
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C.  MATCHED  FILTERS 

Some  interesting  properties  of  matched  filters  will  be 
listed  here.  These  properties  are  derived  in  the  literature 
[Ref.  3,4]. 

A matched  filter  is  the  best  linear  filter  for  detection 
of  a pulse  signal  v(t)  in  noise.  The  impulse  response  h(t) 
of  such  a filter  is  a delayed,  time  inverted  replica  of  the 
input.  If  v(t)  is  the  input  to  the  matched  filter,  its  impulse 
response  is 

h(t)  = Mv[-(t-to)],  t-0 
h(t)  = 0 , t<0. 

where  M = an  arbitrary  constant 


tQ  = time  delay  inherent  in  the  filter. 


From  linear  system  theory,  the  output  g(t)  of  a matched 
filter  is 


g(t) 


mX 


+00 


v(t-x)v  [-(x-t  )]dx 
o 


Let  A = t-x  to  obtain 

* +OO 


gCt) 


yL 


v(\ )v(-t+\+to)d\ 


^ -f  CO 

or  g(t)  = M / v(X)v(X-OdX  = M R CO 

sj  —00 

where  £ = t-to. 

The  output  is  maximized  when  £ = 0 or  t = tQ . 

It  is  concluded  then,  that  the  output  g(t)  of  the  matched 
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filter  is  the  autocorrelation  function  of  the  input.  Fig.  3 
illustrates  the  concept. 


Matched 

F il  ter1 

v(t)' 

i 

h(t) 

i 

Fig.  3.  Matched  filter. 


Matched  filters  for  two-level  voltages  (binary  sequences) 
can  be  realized  by  tapped  delay  lines  or  shift  registers  as 
shown  in  Fig.  4.  The  tapped  delay  line  realization  uses  in- 
verters at  these  positions  where  a zero  element  occurs  in  the 
sequence.  The  shift  register  realization  uses  a reference 
register  where  the  original  sequence  is  stored.  This  can  be 
a read  only  memory  (ROM)  for  example.  Another  register 
receives  the  input  sequence  v(t) . 

In  both  realizations,  +1  units  of  current  flow  through 
the  load  resistor  R,  for  each  element  of  the  input  sequence 
that  agrees  with  the  "stored"  sequence.  And  -1  unit  of  current 
flows  through  R^  for  each  element  of  the  input  sequence  that 
disagrees  with  the  "stored"  sequence.  The  net  output  current 
through  R^  (and  voltage  across  R^)  is  proportional  to  the 
number  of  elements  which  agree  less  the  number  of  elements 
which  disagree.  The  output  g(t)  is,  therefore,  a measure  of 
the  autocorrelation  function  of  the  input  sequence. 

Since  the  systems  of  Fig.  4 perform  a discrete  comparison 
and  siamming  instead  of  multiplication  and  integration,  then 
g(t)  is  a discrete  version  of  Rv(t).  For  example,  the  sequence 
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Fig.  4.  Matched  filter  realizations. 
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1101  has  the  autocorrelation  function  shown  in  Fig.  2,  whereas 
the  output  gd(t)  of  the  discrete  matched  filter  corresponding 
to  1101  is  as  shown  in  Fig.  5. 


tj 

R d 

V 

(t) 

1 

— 4 £ — 3£  —2c 

-0 

£ 

2s  3£  4S 

-i 

Fig.  5.  Discrete  matched  filter  output. 

For  a sequence  of  n elements,  the  peak  output  is  n units 
of  voltage.  It  is  clear  then,  how  signal  detectability 
improves  as  n increases . 

All  binary  sequences  have  autocorrelation  functions  with 
sidelobes  of  various  values.  The  sequences  of  interest  are 
the  ones  with  either  small  or  negative  sidelobes . The 
problem  is  to  find  these  sequences. 

The  next  section  addresses  the  issue  of  forming  a pair  of 
sequences  or  codes,  which  when  processed  with  matched  filters 
and  the  outputs  added,  yield  a waveform  with  one  mainlobe  and 
no  sidelobes.  This  scheme  provides  good  detectability  of 
binary  sequences . 
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III.  COMPLEMENTARY  SEQUENCES 


This  section  reviews  complementary  sequences  and  their 
basic  properties. 

A set  of  complementary  series  is  defined  as  a pair  of 
equally  long,  finite  binary  sequences  which  have  the  property 
that  the  number  of  pairs  of  like  elements  with  any  given 
separation  in  one  series  is  equal  to  the  number  of  pairs  of 
unlike  elements  with  the  same  separation  in  the  other  series . 

For  example  the  two  series: 


A = 00010010 
B = 00011101 

are  complementary.  In  A there  are  three  like  elements  (denoted 
by  l below)  separated  by  one  element. 


0*0*0  1 0*0  1 0 


In  B there  are: 


o o V i W 

three  unlike  elements  (denoted  by  a below)  separated  by  one 
element . 

Similarly  for  all  possible  separations  the  number  of  like 
elements  in  A and  unlike  elements  in  B are  as  follows: 


Separation 

1 

2 

3 

4 

5 

6 

7 

Series  having  the  complementary  property  were  conceived  by 
Marcel  J.  E.  Golay  in  connection  with  the  optical  problem  of 
infrared  multislit  spectrometry  [Ref.  1,5]. 

Complementary  series  have  interesting  autocorrelation 
functions.  If  the  autocorrelation  of  each  sequence  is  taken 
and  these  two  autocorrelations  summed,  the  result  is  zero  for 
all  t except  x = 0 . At  x = 0,  the  sum  is  twice  that  of  either 
sequence.  Therefore,  the  sum  of  the  autocorrelation  functions 
has  one  main  lobe  and  no  sidelobes. 

For  example,  consider  the  sequences  or  codes 

A = 00010010 
3 = 00011101 

Sequence  A has  autocorrelation  function: 

Ra(x)  = 8,  -1,  0,  3,  0,  1,  0,  1,  0. 

Sequence  B has  autocorrelation  function: 

Rb(x)  = 8,  1,  0,  -3,  0,  -1,  0,  -1,  0. 


Number  of 
Likes  in  A 

3 

3 

4 
2 
2 
1 
1 


Number  of 
Unlikes  in  B 

3 

3 

4 
2 
2 
1 
1 
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Their  sum  is: 


E = Ra(t)  + R0(t)  = 16,  0,  0,  0,  0,  0,  0,  0,  0. 


This  property  can  be  treated  in  equation  form  as  follows: 

Let  a^  and  (i  = 1,  2,  , n)  be  the  elements  of 

two  sequences  A and  B each  of  length  n.  Assume  a^  and  b,.  can 
be  either  +1  or  -1.  Then  the  respective  values  of  the  auto- 
correlation functions  will  be 


n-j 

r» 


E a. 

1 1 

ai  + j 

for 

jSO  . 

n-j 

E b. 
i 

X 

bi+j 

for 

j>0. 

Also  cj 

= c • for 

-3 

j<  o . 

d. 

= d . for 

j<0. 

3 

-3 

The  necessary  and  sufficient  condition  for  the  pair  of 
sequences  to  be  complementary  is : 


c . + d . 
3 3 

and  Cj  + dj 


0 for  j / 0 
2n  for  j = 0 . 


When  the  elements  of  the  sequence  are  0 or  1,  then  the 
autocorrelation  function  is  obtained  by  modulo-two  addition. 

In  this  case  the  necessary  and  sufficient  condition  for  the 
series  to  be  complementary  is  that 

n-j  n-j 

I (a. ©a.,.)  = E (b.  © b . ©1)  for  all  j , 1<  j <:  n-1 

i i+j  l i+ j J J 

i=l  i=l 

This  complementary  property  can  be  tested  experimentally 
by  using  matched  filters  since  the  output  of  a filter  matched 
to  its  input  is  the  autocorrelation  function  of  that  input. 
This  realization  is  shown  in  Figure  6 . The  output  of  the 
system  of  Figure  6 is  shown  in  Figure  7 for  the  complementary 
sequences 

A = 111-1 
B = 11-11 

A.  GENERAL  PROPERTIES 


1 . Number  of  Elements 

The  number  of  elements  in  two  complementary  series 
are  equal.  If  they  were  different,  the  pair  of  extreme  elements 
of  the  longest  series  would  remain  unmatched  by  an  unlike  pair 
of  elements  with  the  same  spacing  in  the  other  series . 


2 . Symmetry 

Two  complementary  series (A, B)  are  interchangeable  (3, A); 
that  is,  one  can  take  the  place  of  the  other.  This  results  from 
the  symmetry  of  the  definition  with  respect  to  two  complementary 
series . 
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Fig.  7.  Analog  matched  filter  processing  of 
two  complementary  sequences. 
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Sequence  Length  Even 


A necessary  condition  for  a sequence  pair  to  be  com- 
plementary is  that  their  length  n be  an  even  number. 

4 . Sequence  Length  Sum  of  Two  Squares 

Another  necessary  condition  for  a pair  of  sequences 
to  be  complementary  is  that  their  length  be  the  sum  of  the 
squares  of  two  integers.  The  proof  was  developed  by 
S.  Jauregui  [Ref.  4]: 

5 . Transformations 

A single  pair  of  complementary  series  can  be  the  basis 
for  the  construction  of  64  pairs  of  complementary  series. 

a.  Order  of  Complementary  Sequences 

Denote  the  reverse  of  A by  A . For  example, 
if  A = 1110,  then  A^  = 0111.  The  order  of  the  elements  of 
either  or  both  of  a pair  of  complementary  series  may  be 
reversed.  This  follows  from  the  fact  that  by  reversing  a 
sequence  its  autocorrelation  function  remains  the  same.  The 
proof  is  developed  in  Section  V. 

b.  Complementing  the  Sequence 

Denote  the  complement  of  A by  A.  For  example,  if 
A = 1110,  then  A = 0001.  One  or  both  of  a pair  of  complementary 
sequences  can  be  complemented — putting  zeros  in  the  place  of 
ones  and  ones  in  the  place  of  zeros,  without  affecting  their 
complementary  property.  This  follows  from  the  fact  that  by 


a 
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complementing  a sequence,  its  autocorrelation  function  remains 
the  same.  Section  V provides  the  proof. 

c.  Complementing  Elements  of  Even  Order 

Denote  the  complement  of  the  even  order  elements 
of  A by  A . For  example,  if  A = 1110,  then  Ag  = 1011.  Com- 
plementing the  elements  of  even  order  in  each  sequence — putting 
zeros  in  the  place  of  ones  and  ones  in  the  place  of  zeros,  does 
not  affect  their  complementary  property. 

It  is  concluded  from  the  above  properties  that  a single 
pair  (A, 3)  of  complementary  sequences  can  be  the  basis  for  the 

g 

construction  of  2 =64  pairs  of  complementary  series  (some  of 

which  might  be  identical)  by  either  performing  or  not  performing 
the  following  six  operations: 

a.  Interchanging  the  sequences. 

b.  Reversing  the  first  sequence. 

c.  Reversing  the  second  sequence. 

d.  Complementing  the  first  sequence. 

e.  Complementing  the  second  sequence. 

f.  Complementing  the  elements  of  even  order  of  each 
sequence . 

As  an  example,  consider  the  complementary  pair  A = 00010010 
and  B = 00011101: 

Applying  a gives  B = 00011101  and  A = 00010010. 

Applying  b gives  A^  = 01001000  and  B = 00011101. 

Applying  c gives  A = 00010010  and  B^  = 10111000. 

Applying  d gives  A = 11101101  and  B = 00011101. 
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Applying  e gives  A = 00010010  and  B = 11100010. 

Applying  f gives  Ag  = 01000111  and  Bg  = 01001000. 

By  applying  the  above  properties  properly,  the  original 
pair  can  be  reproduced: 

B = 00011101  and  A = 00010010 

B = 10111000  and  A = 00010010 
r 

3 = 10111000  and  A„  = 01001000 

r r 

B = 01000111  and  A„  = 01001000 
r r 

(f  ) = A = 00010010  and  (A  ) = B = 00011101. 

re  re 

The  last  pair  is  the  same  as  the  original  one . 

6 . Allowable  Lengths 

Since,  as  was  mentioned  before,  the  number  of  elements 
in  complementary  sequences  must  be  even  and  equal  to  the  sum 
of  two  squares,  the  allowable  sequence  lengths  up  to  SO  are 

2,  4,  3,  10,  IS,  13,  20,  26,  32,  34,  36,  40,  SO. 

It  has  been  verified  by  trial,  though,  that  complementary 
sequences  for  length  18  do  not  exist. 

7 . Hamming  Distance 

The  Hamming  distance  of  two  binary  sequences  A and  B 
is  defined  as  the  number  of  positions  in  which  these  two 
binary  sequences  differ.  This  can  be  written  in  modulo  two 
notation  as  follows: 
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i=n 

D( A, B)  = E a.  © b . 

i=l  1 1 

For  example,  the  Hamming  distance  of  the  two  sequences 

A = 0100 
B = 1111 

is 

i = 4 

D(A, B)  = Z a.  ©b.  =(l+0+l+l)=3 

i = l 1 1 

Now  for  a complementary  pair  of  sequences,  it  has  been 
proven  that  their  Hamming  distance  is  always  = j [Ref.  2]. 

For  example,  the  complementary  pair  of  length  n = 10 

A = 1001010001 
B = 1000000110 

has  Hamming  distance 

D(A, B)  =|=5. 

8 . Kernels 

A Kernel  is  a basic  length  sequence  which  cannot  be 
decomposed  into  shorter  length  sequences . The  shortest 
possible  complementary  pair  is  11  and  10.  This  pair  or  any 
of  each  transformation,  which  was  mentioned  before,  is  called 
a kernel  of  length  two  or  a quad. 
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Some  possible  Kernel  lengths  are 


2,  10,  18,  26,  34,  50. 

It  might  be  the  case,  though,  that  complementary  pairs  for 
some  of  them  do  not  exist.  For  example,  lengths  n = 4 and 
n = 8 have  complementary  pairs,  but  are  not  Kernels  because 
they  can  be  constructed  from  n = 2 and  n = 4 sequences 
respectively . 

Among  all  the  above  mentioned  Kernel  lengths,  it  has 
been  verified  by  M.  J.  Golay  that  n = 18  does  not  exist. 

Also,  it  has  been  verified  by  S.  Jauregui  [Ref.  2] 
through  exhaustive  computer  search  that  for  n = 26  only  the 
Kernel  shown  in  Table  I exists,  not  taking  into  account  all 
allowable  transformations. 

For  n = 34  a non-exhaustive  computer  search  by 
S.  Jauregui  revealed  no  Kernel.  An  exhaustive  computer  search 
was  not  possible,  due  to  the  great  computer  time  required. 

It  is  possible,  however,  that  a complete  search  for 
n = 34  could  be  achieved  in  the  future,  using  new  techniques. 
The  following  table  shows  the  Kernels  of  n = 2,  10,  26, 
ignoring  allowable  transformations. 


Table  I 

Kernels  of  Length  2,  10,  26 
Number  of 


n 

Kernels 

A Sequence 

B Sequence 

2 

1 

10 

11 

10 

2 

1001010001 

0101000011 

1000000110 

0000100110 

26 

1 

010011011110101111 

00111010 

101100100001111111 

00111010 
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Note:  (a)  The  possible  Kernel  of  n = 18  does  not  exist. 


(b)  Partial  computer  search  for  n = 34  found  no 
Kernel . 

9 . Number  of  Ones  in  Complementary  Sequences 

S.  Jauregui  showed  [Ref.  2]  that  the  equation 
n = (n-p-q)2  + (p~q)2 

holds  for  two  complementary  sequences  A and  3 of  length  n 
where  p is  the  number  of  ones  in  A and  q the  number  of  ones 
in  B . 

This  leads  to  the  conclusion  that  the  number  of  ones 
in  each  of  the  sequences  of  a complementary  pair  cannot  be 
arbitrary,  but  has  to  satisfy  the  above  relation. 

For  example,  fov  complementary  sequences  of  length 
n = 2 the  number  of  ones  in  A and  B can  be  respectively 


either  (2,  1) 
or  ( 1 , 0 ) . 

3.  SYNTHESIS 

If  the  sequences  A, 3 are  complementary,  they  can  be  used 
to  generate  other  complementary  pairs  as  follows: 

(a)  If  A = aL  a2  a3  — -a^  an 

B = bi  b2  b3  bn_1  bn 

are  a complementary  sequence  pair,  then  the  sequences 


J 
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C = a,  a0 b,  b0 b 

12  n 1 2 n 

D = a,  a0 S.  S~ b„ 

12  n 1 2 n 


are  also  complementary . 

For  example,  consider  the  complementary  pair: 

A = 0001  with  autocorrelation  R^(t)  = 1+5  1*  -1»  0 

B = 0010  with  autocorrelation  Rg(r)  = 4,  -1,  0,  1,  0. 

Then  the  sequences 


C = 00010010 
D = 00011101 


are  also  complementary  with  autocorrelation  functions 
Rc(t)  = 8,  -1,  0,  3,  0,  1,  0,  1,  0 
Rd(t)  = 3,  1,  0,  -3,  0,  -1,  0,  -1,  0. 


(b)  If  A = a.  a0  a,  a , a 

123  n-1  n 


B = b,  b0  b, b_  -i  b 

123  n-l  n 


are  a complementary  sequence  pair,  then  the  sequences 

C = al  b]_  a2  b2 an  bn 

D = al  B1  a2  B2  an  Bn 

are  also  complementary. 

For  example,  consider  the  same  sequences 
A = 0001 
B = 0010 
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The  sequences 


C = 00000110 
D = 01010011 

are  also  complementary  with  autocorrelation  functions 

Rc(t)  = 8,  3,  0,  1,  0,  -1,  0,  1,  0 
Rd(t)  = 8,  -3,  0,  -1,  0,  1,  0,  -1,  0. 

(c)  If  (A,B)  (C,D)  are  two  complementary  sequences 
pairs,  A of  length  n and  C of  length  n,  then  the  pair 

V1  = AC1  AC2 ACn  3dl  Bd2 Bdn 

V2  = Adn Adl  Bcn 3C1 

is  also  complementary,  where  if  an  exponent  is  one  the  A or  B 
sequence  is  left  unchanged,  whereas  if  the  exponent  is  zero 
the  A or  B sequence  is  complemented. 

For  example,  consider  the  complementary  pairs 


A = 11 
B = 10 


and 


C = 00 
D = 01. 


Then  the  pair 

v1  = 00000110 
v 2 = 11001010 
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is  also  complementary  with  autocorrelation  functions 


(t  ) 


(t ) 


8 , 3,  0 , 1,  0 , -1,  0 , 1,  0 . 

8,  -3,  0,  -1,  0,  1,  0,  -1,  0. 


The  above  methods  make  possible  the  generation  of  comple- 
mentary sequence  pairs  of  greater  lengths  than  the  original 
ones.  They  can  be  applied  in  succession  to  generate  very  long 
sequences  which  are  very  useful  in  many  applications. 

For  example  such  a complementary  pair  used  in  a communica- 
tions system  with  matched  filter  processing  like  that  of  Fig.  6 
can  improve  signal  detectability  in  the  presence  of  considerable 
noise  since  the  summer  output  voltage  will  consist  of  a large 
main  lobe  and  no  sidelobes. 


C.  SUPPLEMENTARY  AND  CYCLIC  COMPLEMENTARY  SEQUENCES 

Complementary  sequences  are  subsets  of  two  larger  sets, 
namely  supplementary  and  cyclic  complementary  sequences. 

1 . Supplementary  Sequences 


Consider  two  sequences 


A = al  a2  a3  an_x  an 


B = bl  b2  b3  - -bn_1  bn 
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Now  let 


I = al  a3  a5 an-l 

II  = an  an-2  — “a2 

III  = b,  b3  b5  — -bn_. 


IV 


b 


n 


b 

n- 


2 


The  expression  of  the  sequence  pair  A,  B in  the  form 
(I,  II,  III,  IV)  is  called  sequence  quadruple.  Supplementary 
sequences  are  quadruples  of  sequences  with  the  property  that 
the  total  number  of  likes  at  each  spacing  equals  the  total 
number  of  unlikes  at  the  same  spacing. 

In  terms  of  their  autocorrelation  function,  the  sum 
of  the  four  autocorrelation  functions  is  zero  any  place  but 
t = 0,  where  it  is  four  times  the  length  of  the  sequences. 

For  example, 

A = 1001010001 

B = 1000000110 

Writing  in  (I,  II,  III,  IV)  form  gives 

I = 10000 

II  = 10110 

III  = 10001 

IV  = 01000. 
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The  new  sequences  (I,  II,  III,  IV)  have  autocorre- 
lation functions 


Rx(t)  = 5,  2,  1,  0,  -1,  0 
Ri;[(t)  = 5,  -2,  -1,  2,  -1,  0 
Riix(t)  = 5,  0,  -1,  -2,  1,  0 
RIV(x)  = 5,  0,  1,  0,  1,  0. 

The  sum  of  the  autocorrelation  functions  is 

I = Rj(t)  + Rxt(t)  + + RjyCx) 

= 20,  0,  0,  0,  0,  0. 

In  this  example  the  sequences  A,  3 are  complementary. 
In  general  they  do  not  have  to  be  though,  since  the  supplementary 
sequences  are  a larger  set.  This  is  illustrated  in  the  following 
example . 

Consider 

I = 000100111 
II  = 000101001 
III  = 000101000 
IV  = 000110110 

with  autocorrelation  functions 

Rx(t)  = 9,  2,  -1,  0,  1,  0,  -3,  -2,  -1,  0 

RI];(t)  = 9,  -2,  1,  0,  -1,  2,  1,  0,  -1,  0 

RIII(t)  = 9’  3’  ~2’  1’  °»  3’  2’  0 

Riv(t)  = 9,  0,  -3,  2,  -1,  -2,  -1,  0,  1,  0. 
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The  sum  is 


Z = R.j.(t)  + Rjj(t)  + RII];(t  ) + Riv(t) 
= 36,  0,  0,  0,  0,  0,  0,  0,  0,  0 


which  is  zero  except  at  the  position  t = 0,  where  it  is  36. 

From  the  sequences  (I,  II,  III,  IV)  the  sequences 
A,  B can  be  constructed 


A = 010000110001101010 


3 = 000101100111000000 


In  this  example  neither  the  A,  B sequences  nor 
the  (I,  II,  III,  IV)  ones  are  complementary,  which  demonstrates 
the  fact  that  supplementary  sequences  are  a larger  set. 


Figure  8 gives  another  example  of  the  supplementary 


property . 


2 . Cyclic  Complementary  Sequences 


In  general  a cyclic  sequence  is  a never  ending  periodic 
sequence  of  zeros  and  ones  which  has  period  of  n elements.  A 
cyclic  complementary  sequence  pair  is  a pair  of  cyclic 
sequences,  each  of  period  n,  where  the  number  of  likes  in  one 
sequence  equals  the  number  of  unlikes  in  the  other  one,  for 
all  possible  spacings . In  equation  form 

Ci  ®ViH  bl  ® bn-i +1  ffil- 


In  terms  of  autocorrelation  functions  the  sum  of  the 
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two  periodic  autocorrelation  functions  is  zero  except  at 
T = kn,  where  k = 0,  1,  2,  .... 


For  example,  the  cyclic  complementary  sequences 

A = 01010011,  01010011,  

B = 00000110,  00000110,  

have  periodic  autocorrelation  functions 

Ra(t)  = 8,  -4,  0,  0,  0,  0,  0,  -4,  8,  

Rg(t ) = 8,  4,  0,  0,  0,  0,  0,  4,  8,  .... 

The  sum  is 

Z = Ra(t)  + RgCx)  = 16,  0,  0,  0,  0,  0,  0,  0,  16 

Figure  9 demonstrates  the  above  example. 

Complementary  sequences  are  always  supplementary  and 
cyclic  complementary,  but  the  opposite  is  not  always  true. 
Supplementary  and  cyclic  complementary  sequences  constitute 
a larger  set. 
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Fig.  9 


. The  cyclic  complementary  property. 
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Complementary  sequences  are  attractive  for  use  in  communi- 
cations and  ranging  systems  because  the  usable  receiver  output 
voltage  has  no  sidelobes.  Sequences  providing  "small"  side- 
lobes  compared  to  the  main  lobe  may  also  be  useful  in  some 
applications . 

It  may  even  be  desirable  to  have  a small  sidelobe  at  a 
known  position  and  level.  For  example,  such  a sidelobe  can 
be  used  to  measure  doppler  as  shown  in  Section  VII. 

Two  binary  sequences  whose  summed  autocorrelation  functions 
exhibit  two  small  sidelobes  are  called  almost  complementary  in 
this  report,  since  they-  exhibit  properties  similar  to  comple- 
mentary sequences. 

It  should  be  noted  that  there  are  many  sequences  with  small 
sidelobes,  but  the  ones  of  interest  here,  are  those  for  which 
certain  rules  hold.  By  applying  these  rules,  the  magnitude 
and  position  as  well  as  polarity  (positive  or  negative)  of  the 
sidelobe  can  be  precisely  predicted.  By  knowing  these  rules 
almost  complementary  sequences  or  codes  can  be  constructed. 

An  extensive  search  for  "good"  almost  complementary 
sequences  was  made.  As  a result,  it  is  possible  to  list  general 
rules  for  constructing  almost  complementary  sequences  with  pre- 
dictable sidelobe  levels,  polarity  and  position.  There  is  a 
distinction  between  positive  sidelobe  and  negative  sidelobe 
sequences . 


J 
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A.  POSITIVE  SIDELOBE  SEQUENCES 


Here  two  cases  are  developed. 
+ 1 

1.  Sidelobes  - j K Away 


If  A,  B are  complementary  sequences  of  length  n and 
the  new  sequences  C and  D of  length  K = 2n  are  formed  as 
follows , 

C = AA 
D = BB 


the  sum  of  the  autocorrelation  functions  of  C and  D gives  only 
two  positive  sidelobes  of  magnitude  half  that  of  the  main  lobe 
at  t = - j K. 

For  example,  if  n = 8 , then  K = 16  with  the  main  lobe 
level  of  16,  and  sidelobe  levels  of  8 at  t = - 4, 

For  example,  consider  the  two  complementary  sequences 
of  length  n = 4 , 

A = 0001 
B = 1011, 

and  construct  the  almost  complementary  sequences 

C = AA  = 00010001 
D = BB  = 10111011 

having  autocorrelation  functions 


L 


Rc(t)  = 8,  1,  0,  -1,  4,  1,  0,  -1,  0 
Rd(t)  = 8,  -1,  0,  1,  4,  -1,  0,  1,  0. 

Adding  gives 

Z = Rc(t)  + Rd(t)  = 16,  0,  0,  0,  8,  0,  0,  0,  0. 

So  there  are  only  two  positive  sidelobes  with  levels  half  the 

+ 1 + 

level  of  the  main  lobe  and  at  a distance  t = - K = - 4 . 

Figure  10  illustrates  the  above  example. 

Another  example  is 

A = 0001001011100010 
B = 0100011110110111 

A,  B are  complementary  of  length  n = 16.  Constructing  C = AA, 

D = BB  gives 

C = 00010010111000100001001011100010 
B = 01000111101101110100011110110111 

C,  B are  almost  complementary  of  length  K = 32  with  auto- 
correlations 

Rc(t)  = 32,  -1,  0,  -1,  0,  3,  0,  -5,  0,  -1,  0,  9,  0, 
1,  0,  1,  16,  -1,  0,  -1,  0,  -1,  0,  -1,  0,  -3, 
0,  5,  0,  1,  0,  1,  0 . 

Rd(t)  = 32,  1,  0,  1,  0,  -3,  0,  5,  0,  1,  0,  -9,  0, 

-1,  0,-1,  16,  +1,  0,  +1,  0,  +1,  0,  +1,  0,  3, 

0 , -5 , 0 , -1 , 0 , -1 , 0 
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and  sum 


E = Rc(t)  + RdCt)  = 64,  0,  0,  0,  0,  0,  0,  0,  0,  0,0, 

0,  0,  0,  0,  0,  32,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0,  0 

with  main  lobe  level  64  and  positive  sidelobe  levels  32  at 
t = — 16  as  predicted  by  the  rule. 

For  this  example,  the  IBM- 360  computer  was  used. 

2 

2 . Sidelobes  - y K Away 

Let  A,  B be  complementary  sequences  of  length  n.  Let 

| be  a new  sequence  generated  by  taking  only  the  j first  digits 

of  A (truncating  A after  its  y first  digits).  Also  let  $ be  a 

new  sequence  generated  by  taking  only  the  y first  digits  of  B 

(truncating  B after  its  y first  digits).  Also  let  $ be 

complementary  of  length  y.  Then  the  new  sequences  of  length 
n 3n 

K = n + y - ~2~  are 

C = A^ 

D = B$ 

and  C and  D are  almost  complementary  sequences.  The  sum  of 

their  autocorrelation  functions  has  two  positive  sidelobes 

1 + 2 
only  of  magnitude  y that  of  the  main  lobe  and  at  t = - y K. 

In  this  case  there  is  an  improvement  relative  to  the  previous 

case  in  that  the  sidelobes  are  smaller  compared  to  the  main 

lobe  and  farther  removed  from  the  main  lobe . 
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The  sequences  A,  B of  length  n can  be  constructed  by  using 
two  complementary  sequences  of  length  y according  to  the  rule 
in  Section  III.B(a).  Then  the  sequences  $ are  always 
complementary.  For  example,  consider  the  two  complementary 
sequences  of  length  n = 8 

A = 00011101 
B = 01001000. 

The  new  sequences  of  length  12  are 

C = A|  = 000111010001 
D = B$  = 010010000100 

with  autocorrelation  functions 

Rc(t)  = 12,  1,  0,  -5,  0,  -5,  0,  1,  4,  1,  0,  -1,  0. 
Rd(t)  = 12,  -1,  0,  5,  0,  5,  0,  -1,  4,  -1,  0,  1,  0. 

and  sum 

Z = Rc(t)  + Rd(t)  = 24,  0,  0,  0,  0,  0,  0,  0,  8,  0, 

0,  0,  0 

with  main  lobe  level  24  and  sidelobe  levels  "3  (24)  = 8 at  a 
distance  r = - 8 . 

This  example  is  illustrated  in  Fig.  11.  Another 
example  is 

A = 01001000000111010100100011100010 
B = 01001000000111011011011100011101 
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Fig.  10.  Positive  sidelobes  at  t = - j K. 
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A,  B are  complementary  of  length  n = 32. 

The  new  sequences  of  length  K = 48  are 

C = A A = 01001000000111010100100011100010 
0100100000011101 

D = B$  = 1011100011101101101x100000010010 
1011100011101101 

with  autocorrelation  functions 

Rc(t)  = 48,  1,  2,  3,  0,  -1,  -2,  1,  0,  3,  -2,  5,  0, 
+9,  2,  -5,  0,  -5,  2,  9,  0,  5,  -2,  3,  0,  1, 

-2,  -1,  0,  3,  2,  1,  16,  1,  2,  3,  0,  -1,  -2, 

1,  0,  -1,  -2,  1,  0,  -3,  2,  -1,  0. 

R^Ct)  = 48,  -1,  -2,  -3,  0,  1,  2,  -1,  0,  -3,  2,  -5, 

0,  -9,  -2,  5,  0,  5,  -2,  -9,  0,  -5,  2,  -3, 

0,  -1,  2,  1,  0,  -3,  -2,  -1,  16,  -1,  -2,  -3, 
0,  1,  2,  -1,  0,  1,  2,  -1,  0,  3,  -2,  1,  0. 

and  sum 

Z = Rc(t)  + Rd(t)  = 96,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

32,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0. 

The  main  lobe  has  level  96  and  the  sidelobes  have  levels 

2 

at  a distance  r = - j 48  = - 32  as  predicted  by  the  rule. 
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rr I 

Here  also  the  IBM- 360  computer  was  used  because  of  the 
length  of  the  sequence. 

[ 

B.  NEGATIVE  SIDELOBE  SEQUENCES 

Here  two  cases  are  developed  similar  to  the  ones  considered 
before . 

+ 1 

1.  Sidelobes  - j K Away 


If  A,  B are  complementary  sequences  of  length  n and  the 
new  sequences  C = AA  and  D = BB  of  length  K = 2n  are  formed, 
where  A,  B represent  the  complements  of  A,  B respectively,  then 
the  sum  of  their  autocorrelation  functions  gives  two  negative 
sidelobes  only  with  magnitudes  half  that  of  the  main  lobe  and 
at  a distance  x = - ^ K. 

For  example 

A = 1011 
B = 1110 

are  complementary.  Then 

C = AA  = 10110100  and 
D = BB  = 11100001 

are  almost  complementary  with  autocorrelation  functions 


Rc(x)  = 8,  -3,  0,  3,  -4,  1,  0,  -1,  0 
Rd(t)  = 8,  3,  0,  -3,  -4,  -1,  0,  1,  0 
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and  sum 


Z = Rc(t)  + Rd(t)  = 16,  0,  0,  0,  -8,  0,  0,  0,  0. 

Here  the  sidelobes  are  negative  with  levels  half  that  of  the 
main  lobe  at  t = - | K = - j 9 = - 4.  This  is  illustrated 
in  Fig.  12. 

Another  example,  for  which  the  IBM- 360  computer  was 

used  is: 


A = 0001001011100010 
B = 0100011110110111 


A,  B are  complementary  of  length  n = 16 . Constructing 
C = AS,  D = B5  gives 

C = 00010010111000101110110100011101 
D = 01000111101101111011100001001000 


Sequences  C and  D are  almost  complementary  of  length 
K = 32  with  autocorrelation  functions 


Rc(t)  = 32,  -3,  0,  -3,  0,  -7,  0,  1,  0,  -5,  0,  11,  0, 
3,  0,  3,  -16,  1,  0,  1,  0,  1,  0,  1,  0,  3,  0, 
-5  , 0 , -1 , 0 , -1 , 0 


Rd(t) 


32,  3,  0,  3,  0,  7,  0,  -1,  0,  5,  0,  -11, 
0,  -3,  0,  -3,  -16,  -1,  0,  -1,  0,  -1,  0, 


-1,  0,  -3,  0,  5,  0,  1,  0,  1,  0 


\ 
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and  sum 


L 


Z = Rc(t)  + Rd(t)  - 64,  0,  0,  0,  0,  0,  0,  0,  0,  0, 


0, 

0, 

0, 

0, 

0, 

0,  -32,  0 

0, 

0, 

0, 

0, 

0, 

0,  0,  0, 

0,  0 


with  main  lobe  level  64  and  negative  sidelobe  levels  32  at 
x = - 16  as  predicted  by  the  rule. 

2 

2.  Sidelobes  - y K Away 

Let  A,  B be  complementary  sequences  of  length  n.  Let 
J be  a new  sequence  generated  by  taking  only  the  complement  of 
the  y first  digits  of  A. 

Let  $ be  a new  sequence  generated  by  taking  only  the 
complement  of  the  y first  digits  of  B.  Also,  let  $ be 
complementary  of  length  Then  the  new  sequences  C = A$  and 

3 T"1 

D = Bp  of  length  K = n + '£  = — will  be  almost  complementary 

and  the  sum  of  their  autocorrelation  functions  will  have  two 

negative  sidelobes  only  with  magnitudes  one-third  that  of  the 

. + 2 
main  lobe  at  x = - | K. 

For  example, 

A = 11101101 
B = 10111000 

are  complementary  sequences  of  length  n = 8 . The  new  sequences 
C = A$  = 111011010001  and  D = B|  = 101110000100  of  length  K = 12 
are  almost  complementary  with  autocorrelations 
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Rc(t)  = 12,  -1,  0,  1,  0,  -1,  0,  1,  -4,  -1,  0,  1,  0 

Rd(t)  = 12,  1,  0,  -1,  0,  1,  0,  -1,  -4,  1,  0,  -1,  0 

and  sum 

£ = Rc(t)  + Rd(t)  = 24,  0,  0,  0,  0,  0,  0,  0,  -8,  0, 

0,  0,  0 

The  only  sidelobes  are  negative  with  level  ^ (24)  = 8 at 

t = - j K = - ~ (12)  = - 8.  This  is  illustrated  in  Fig.  13. 

Another  example  obtained  with  the  use  of  the  IBM-360 
computer  is 


A = 01001000000111010100100011100010 
B = 01001000000111011011011100011101 

A,  B are  complementary  of  length  n = 32 . The  new  se- 
quences of  length  K = 48  are 

A = A$  = 0100100000011101010010001110001010110 
11111100010 

B = Bf  = 1011100011101101101110000001001001000 
11100010010 . 

The  sequences  C and  D are  almost  complementary  with  auto- 
correlations 

Rc(t)  = 48,  -1,  6,  -3,  0,  1,  -6,  -1,  0,  1,  -6,  -1, 

0,  3,  6,  1,  0,  -1,  -6,  -3,  0,  1,  6,  -1,  0, 

1,  6,  -1,  0,  3,  -6,  1,  -16,  -1,  -2,  -3,  0, 

1,  2,  -1,  0,  1,  2,  -1,  0,  3,  -2,  1,  0 
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RjjCt)  = 48 , 1,  -6,  3,  0,  -1,  6,  1,  0,  -1,  6,  1,  0, 


-3,  -6,  -1,  0,  1,  6,  3,  0,  -1,  -6,  1,  0, 
-1,  -6,  1,  0,  -3,  6,  -1,  -16,  1,  2,  3,  0, 
-1,  -2,  1,  0,  -1,  -2,  1,  0,  -3,  2,  -1,  0 


and  sum 


I = Rc(t)  + Rd(t)  = 96,  0,  0,  0,  0,  0,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 
0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 
0,  0,  -32,  0,  0,  0,  0,  0,  0, 

0,  0,  0,  0,  0,  0,  0,  0,  0,  0 

with  a main  lobe  level  of  96  and  negative  sidelobe  levels  of 
32  at  r = - 32  as  predicted  by  the  rule. 

In  signal  detection  applications  almost  complementary 
sequences  with  negative  sidelobes  offer  better  noise  immunity 
than  those  with  positive  sidelobes.  The  negative  sidelobes 
can  be  removed  with  an  envelope  detector. 
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V.  NUMBER  OF  DIFFERENT  AUTOCORRELATION  FUNCTIONS  IN  ALL 

BINARY  SEQUENCES  OF  FIXED  LENGTH 

Binary  sequences  useful  for  communications  or  ranging 
purposes  have  autocorrelation  functions  with  small  sidelobes. 
Small  can  be  defined  in  terms  of  a predetermined  level. 

How  can  these  sequences  be  found?  At  present,  the  only 
way  to  find  these  sequences  is  to  form  the  autocorrelation 
functions  of  all  possible  sequences  of  a given  length  and 
then  select  the  desirable  ones . This  is  a tedious  task 
specially  for  long  sequences,  since  there  are  2n  possible 
binary  sequences  of  length  n.  However,  as  shown  in  this 
section,  many  of  these  2n  sequences  have  the  same  auto- 
correlation function. 

In  general  there  are  four  sequences  of  given  length  which 
have  the  same  autocorrelation  function: 

1.  The  sequence  itself. 

2.  The  sequence  obtained  by  reversing  the  original 
sequence . 

3.  The  sequence  obtained  by  complementing  the  original 
sequence . 

4.  The  sequence  obtained  by  complementing  and  reversing 
the  original  sequence. 

For  example,  consider  the  sequence  of  length  n = 4: 

A = 1101 
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Reversing  A gives 


A = 1011  = B. 
r 

If  A is  complemented,  results  in 

A = 0010  = C. 

If  A is  reversed  and  complemented,  gives 

A = 0100  = D. 
r 

All  these  sequences  A,  B,  C,  and  D have  the  same  auto- 
correlation function 

RA(t)  = VT)  = Rc(T)  = RD(t)  = 4’  -1»  °»  1>  °- 

These  results  are  easy  to  prove.  Consider  each  case 
separately . 

1.  When  a sequence  is  reversed  and  its  autocorrela- 
tion function  taken,  this  is  exactly  the  same  as  if  the  auto- 
correlation function  of  the  original  sequence  was  taken, 
since  the  formation  of  the  autocorrelation  function  can  be 
considered  as  being  accomplished  by  "sliding"  the  sequence 
past  itself  either  to  the  right  or  to  the  left.  So,  "sliding" 
to  the  right  for  A is  equivalent  to  "sliding"  to  the  left 

for  Ar. 

2.  The  autocorrelation  function  of  a sequence  is 
generated  by  forming  products  and  adding  them.  The  general 
form  of  one  of  these  products  is 
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a^  and  a^  can  have  the  values 


a.  = 1 or  -1 


= 1 or  -1. 


Here  the  values  1 and  -1  are  used  since  multiplication  is 
considered  in  the  formation  of  the  autocorrelation  function 
(if  modulo  2 addition  were  considered,  the  values  1 and  0 
would  be  used) . 

Now  the  possible  values  of  the  product  a^  a^  are 

ai  a^  = (1)(1)  = 1 
or  ai  a j = (1) (-1)  = -1 
or  a^  a j = (-1) (1)  = -1 
or  ai  a j = (-1)(-1)  = 1. 

If  the  sequence  is  complemented,  1 is  replaced  by  -1,  and 
-1  by  1.  So  the  possible  values  of  the  product  a^  a^  are 
respectively 


a . 

l 

a . 

: 

= (-D(-l)  = 1 

or 

a . 

i 

a j = (1)(-1)  = -1 

or 

a . 

i 

a.  = (1)(-1)  = -1 

or 

a . 

i 

a.  = (1)(1)  = 1. 

The  values  of  the  product  are  the  same  as  before 

for  all  i and  j and  so,  the  autocorrelation  function  remains 
the  same . 
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3.  Since  the  autocorrelation  function  is  the  same  if 
the  sequence  is  reversed  or  complemented,  it  follows  that  it 
will  be  also  the  same  if  the  sequence  is  reversed  and  comple- 
mented . 

Some  sequences  are  their  own  reverse.  These  sequences 
are  called  here  symmetric  and  will  be  denoted  by  the  letter  S. 

For  example  the  sequence  A = 1001,  when  reversed  gives 
Ar  = 1001  = A.  In  this  case  only  the  sequence  itself  and  its 
complement  need  to  be  considered.  The  sequences  A = 1001  and 
A = 0110  have  the  same  autocorrelation  function. 

For  some  sequences,  complementing  and  reversing  gives  the 
same  sequence.  These  sequences  are  called  here  R sequences. 

For  example,  the  sequence  A = 000111  when  reversed  and 
complemented  gives  = 000111  = A.  In  this  case  only  the 
sequence  itself  and  its  reverse  need  to  be  considered.  The 
sequences  000111  and  111000  have  the  same  autocorrelation 
function . 

It  is  assumed  that  the  number  of  sequences  having  the 
same  autocorrelation  function  is  either  two  or  four.  It  has 
not  been  shown  though  that  there  does  not  exist  any  other 
number  of  sequences  such  as  3 or  5 or  6 , etc.,  that  have  the 
same  autocorrelation  function.  In  this  work,  it  has  been 
verified  by  exhaustive  computer  search  for  sequences  of  various 
lengths,  that  the  number  of  sequences  having  the  same  auto- 
correlation function  is  either  two  or  four. 

When  an  R sequence  is  reversed  and  complemented,  the  ones 
of  the  R sequence  become  the  zeros  of  the  reverse  complement. 
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Sequence 

Autocorrelation 

000 

3, 

2,  1,  0 

001 

3, 

0,  -1,  o 

010 

3, 

-2,  1,  0 

011 

3, 

O 

M 

O 

100 

3, 

0,  -1,  0 

101 

3, 

-2,  1,  0 

110 

3, 

0,  -1,  o 

111 

3, 

2,  1,  0 

Here  there  are  four  sequences  which  are  neither  R nor  S, 
001,  Oil,  100,  110  with  the  same  autocorrelation  function 
3,  0,  -1,  0.  Also  there  are  four  S sequences:  000,  111,  010, 
101.  The  sequences  000  and  111  have  the  same  autocorrelation 
function  3,  2,  1,  0 and  the  sequences  010  and  101  have  the 
same  autocorrelation  function  3,  -2,  1,  0.  In  this  example 
there  is  no  R sequence  since  n = 3 is  odd. 

Next  the  exact  number  of  R and  S sequences  will  be 
established  among  all  the  possible  sequences  of  length  n. 

The  maximum  number  of  binary  sequences  of  length  n is  given 
by  2n.  The  sequences  of  n = 1 are  0 and  1.  The  number  of 
S sequences  is  two.  So  S = 2 here;(0,  1).  For  n = 2 all 
sequences  are  00,  01,  10,  11.  Here  S = 2:  (00,  11),  R = 2: 

(01,  10),  and  S + R = 4.  For  n = 3,  the  possible  sequences  are 
000,  001,  010,  011,  100,  101,  110,  111.  Here  S = 4:(000,  010, 
101,  111),  R = 0 and  S + R = 4.  Going  to  n = 4 , it  can  be 
seen  that  S = 4 , R = 4 and  S + R = 8 . 
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The  relation  between  n and  S and  R can  be  derived  by 
considering  the  mechanism  of  moving  from  a sequence  of  length 
n to  the  next  one  of  length  n + 1 . 

For  example,  the  sequences  of  n = 2 are  formed  by  taking 
the  sequences  of  length  n = 1 and  adding  in  front  of  each  of 
them  a zero  and  a one,  one  at  a time,  so  the  number  of  sequences 
for  n = 1 is  doubled  and  all  the  possible  sequences  of  length 
n = 2 are  formed.  By  doing  so,  the  following  can  be  noted: 

Ca)  When  moving  from  n even  to  n+1  which  is  odd,  the 
S + R sequences  in  n is  equal  to  the  S sequences  in  n+1. 

(b)  When  moving  from  n odd  to  n+1  which  is  even,  the 
S in  n is  half  the  S + R in  n+1. 

For  example,  for  n = 4 there  are  S + R = 8 sequences  and 
for  n = 5 there  are  S = 8,  but  moving  to  n = 6 gives  S + R = 16 . 

So,  moving  from  an  even  length  to  the  next  keeps  the  number 
S + R,  but  moving  from  an  odd  length  to  the  next  doubles  the 
S + R number. 

Using  formulas  it  can  be  written: 

(a)  n = even,  (S  + R)  = (S)  ,, 

n n+1 

(b)  n = odd,  (S)n  = j (S  + R)n+1 

Now,  considering  the  following  short  table  giving  the  S + R 
terms  of  all  possible  sequences  up  to  length  n = 10 , the 
relation  between  n and  S + R can  be  obtained. 


functions.  So  finally,  all  the  2n  sequences  have 


n+2  n+2 

2n  - 2 2 + 2 2 


different  correlation  functions . By  re- 


arranging this  result,  the  final  formula  for  n even  is  derived 


(T)  = 2n“2  + 2n/2 

even  j 


(b)  n is  odd.  Similarly  all  possible  sequences  here 


number  2n.  Also,  S + R = 2 2 . Deducting  the  2 2 S + R 


n+1 

2n  - 2 2 


sequences  from  2 gives  2-2  sequences  with  

different  autocorrelation  functions.  So,  all  the  2n  sequences 


n+1  n+1 

9n  9 2 9 2 

have  + in- 


different autocorrelation  functions 


or  in  the  final  form 


2n-2  + — 


2 


Two  examples  are  taken 


(a)  n = 7 


(T)odd 


So  the  2 = 128  possible  sequences  of  length  7 give  only  36 

different  autocorrelation  functions. 


(b)  n = 10 


o o ^ 

2 0 + 4j-  = 272 


VI.  COMPUTER  SEARCH  FOR  "GOOD  CODES 


In  this  section  the  computer  programs  for  obtaining  sequences 
with  "small"  autocorrelation  sidelobes  are  discussed.  These 
sequences  or  "good"  codes  are  found  by  applying  mainly  the 
results  of  the  previous  section.  Because  of  the  large  number 
of  different  sequences  of  length  n for  even  modest  values  of  n , 
it  is  necessary  to  use  a digital  computer  to  search  for  "good" 
codes.  Practically,  an  exhaustive  search  is  limited  to  n s 20 
with  present  digital  computers. 

For  small  lengths  (up  to  n = 10)  it  is  possible  that  a 
programmable  calculator  can  be  used  to  find  the  autocorrelation 
function  of  one  sequence  at  a time.  This  was  done  with  a TI-59 
programmable  hand  calculator  by  storing  each  element  of  the 
sequence  in  a memory  location.  Then  the  autocorrelation  function 
was  formed  by  multiplication  of  the  proper  elements  in  each 
position.  In  this  case,  +1  and  -1  is  used  for  the  elements 
of  the  sequence. 

The  algorithm  for  computing  the  autocorrelation  function 
of  a sequence  on  a large  computer  (IBM-360)  is  constructed. 

This  is  given  in  Program  1 on  page  111.  In  this  program  the 
autocorrelation  function  of  only  one  sequence  can  be  computed. 

The  sequence  has  to  be  punched  on  a computer  card.  A small 
modification  gives  Program  2 on  page  112  which  gives 
the  option  to  find  the  autocorrelation  functions  of  any  number 
of  sequences.  Each  of  these  sequences  has  to  be  punched  on 
a separate  card. 
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In  these  computer  programs  and  also  all  the  next  ones,  the 
sequences  are  represented  with  zeros  and  ones . The  algorithm 
for  finding  the  autocorrelation  function  has  been  made  by 
comparing  the  number  of  like  and  unlike  elements  in  every 
position . 

To  avoid  punching  the  sequences  on  the  cards,  a program  is 
created  to  generate  automatically  all  the  possible  2n  sequences 
at  length  n.  This  is  accomplished  by  counting  in  binary  from 
0 to  2n  and  thus  generating  all  the  binary  numbers  of  length  n. 
This  program  is  combined  with  the  program  for  the  computation 
of  the  autocorrelation  function.  So,  every  time  a sequence  is 
generated,  its  autocorrelation  function  is  formed.  This  is 
Program  3 on  page  113. 

Since  only  the  different  autocorrelation  functions  are  of 
interest  here,  a program  is  written  to  select  only  those  codes 
having  different  autocorrelation  functions.  To  understand  the 
operation  of  this  program,  consider  the  following  example. 

Take  all  the  sequences  of  length  n = 4 and  the  corresponding 
autocorrelation  functions . 
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Autocorrelation 


Sequence 


0000 

4, 

3,  2,  1,  0 

0001 

4, 

1,  o,  -1,  o 

0010 

4, 

-1,  0,  1,  0 

0011 

4, 

1,  -2,  -1,  0 

0100 

4, 

-1,  o,  1,  0 

0101 

4, 

-3,  2,  -1,  0 

0110 

4, 

-1,  -2,  1,  0 

0111 

4, 

1,  o,  -1,  0 

1000 

4, 

1,  0,  -1,  0 

1001 

4, 

-1,  -2,  1,  0 

1010 

4, 

-3,  2,  -1,  0 

1011 

4, 

-1,  0,  1,  0 

1100 

4, 

1,  -2,  -1,  0 

1101 

4, 

-1,  0,  1,  0 

1110 

4, 

1,  0,  -1,  0 

1111 

4, 

3,  2,  1,  0 

4 

There  are  a total  of  2 =16  sequences . After  the  eighth 

sequence  0111,  the  other  sequences  are  complements  of  the 
first  eight  ones.  So,  they  give  no  new  autocorrelation 
function,  and  therefore  they  do  not  need  to  be  taken  under 
consideration.  In  the  first  eight  sequences  there  are  two 
pairs  with  the  same  autocorrelation  function  0001  and  0111 
and  also  0010  and  0100.  So,  only  six  different  autocorrela- 
tion functions  remain.  The  general  way  to  proceed  is  to  take 

2n 

each  sequence  in  the  first  -j— , reverse  it,  complement  it  and 
reverse  complement  it  and  then  keep  only  the  original  sequence 
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and  reject  the  others,  since  they  have  the  same  autocorrela- 
tion function.  Here  two  cases  are  considered. 

(a)  The  sequence  ends  in  0. 

In  that  case  its  complement  and  reverse  comple- 
ment will  start  with  one,  so  they  belong  to  the  sequences  after 
2n 

the  first  y ones  and  need  not  to  be  considered.  In  that  case 
only  the  reverse  of  the  sequence  is  taken. 

(b)  The  sequence  ends  in  1. 

In  that  case  its  reverse  and  complement  will  start 

2n 

with  one,  so  they  belong  to  the  sequences  after  the  first  y ones 
and  need  not  to  be  considered.  In  that  case  only  the  reverse 
complement  is  taken. 

Now,  the  program  takes  each  sequence  as  it  is  generated  and 

reverses  it  if  it  ends  in  0 , or  reverse  complements  if  it  ends 

in  1.  If  the  resulting  sequence  represents  a smaller  binary 

number  than  the  original  one,  this  means  '•'hat  the  resulting 

sequence  was  generated  before  and  its  autocorrelation  function 

already  taken,  so  there  is  no  need  to  be  taken  again,  and  the 

program  goes  to  the  next  sequence.  The  procedure  is  repeated 

2n 

until  all  the  first  y sequences  are  finished.  This  way  only 
the  codes  having  different  autocorrelation  functions  are  listed 
in  the  printed  output. 

For  example,  in  the  case  of  n = 4 the  result  is  as  follows: 
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Sequence 


Autocorrelation 


0000 

3,  2,  1,  0 

0001 

4, 

1,  0,  -1,  0 

0010 

4, 

-1,  0,  1,  0 

0011 

4, 

1,  -2,  -1,  0 

0101 

4, 

-3,  2,  -1,  0 

0110 

4, 

-1,  -2,  1,  0 

And  the  number  of  different  autocorrelation  functions  is  6, 
which  is  in  agreement  with  the  formula 


(A) 


even 


2 


n- 


2 


2n/2 

2 


6 


The  program  which  generates  automatically  the  sequences 
of  length  n and  computes  only  the  different  autocorrelation 
functions  is  Program  4 on  page  114 . An  example  for  n = 10  is 
given  on  page  84. 

All  the  different  autocorrelation  functions  are  not  needed. 
Only  those  with  small  sidelobes.  So  a filtering  procedure  has 
to  be  introduced  in  the  program  to  keep  only  those  auto- 
correlations which  have  sidelobes  equal  or  smaller  than  a pre- 
determined level.  Program  5 on  page  115  generates  automatically 
the  sequences  of  length  n and  prints  only  those  different  auto- 
correlation functions  with  sidelobe  levels  equal  or  less  than  1. 
For  this  case  the  lengths  n = 10,  11,  12,  13,  15,  17  and  20 
were  examined,  and  the  number  of  different  autocorrelation 
functions  found  with  sidelobe  levels  equal  or  less  than  1 were, 
respectively,  11,  1,  16,  31,  0,  40,  3. 


73 


It  is  interesting  to  note  that  for  n=ll  there  is  only  one 
such  autocorrelation  which  is  a Barker  code  (a  sequence  with 
sidelobes  between  +1  and  -1) . For  n=15  no  such  codes  exist, 
and  for  n=20  only  three  were  found. 

For  the  case  n=20  it  was  not  possible  to  make  an  exhaustive 
search  because  of  the  computer  time  required.  Using  30 
minutes  of  computer  time,  only  three  such  autocorrelation 
functions  were  found.  It  is  estimated  that  about  2 hours  of 
computer  time  is  required  to  make  an  exhaustive  search. 

It  is  evident  that  for  lengths  greater  than  n=20  even  with 
a large  computer  an  exhaustive  search  is  impractical. 

There  is  though  a way  to  search  regions  of  big  sequences . 
Computer  Program  6 on  page  116  was  used  to  search  for  a region 
of  the  n=20  case  with  no  new  results. 

The  results  for  the  cases  n=10,  11,  12,  13,  17,  20 
are  listed  on  pages  96,  98,  99,  101,  105,  110,  respectively. 
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VII.  APPLICATIONS 


Complementary  sequences,  almost  complementary  sequences 
and  codes  with  small  sidelobes  can  be  used  in  communications, 
ranging  and  spread  spectrum  systems.  This  section  lists  some 
possible  applications. 

Two  complementary  sequences  can  be  transmitted  simultane- 
ously using  a quadriphase  phase-shift  keying  (QPSK),  for 
example.  Two  complementary  sequences  A and  B are  applied  to 
a QPSK  modulator.  The  output  V (t)  is  a sine  wave  with  a 
phase  which  can  have  one  of  four  values.  In  the  receiving 
system,  Vc(t)  is  demodulated.  The  outputs  of  the  demodulator 
are  the  sequences  A and  B.  Each  sequence  is  processed  by  a 
matched  filter  and  the  outputs  of  the  two  matched  filters 
added.  Then  the  output  of  the  summer  will  have  one  main  lobe 
and  no  sidelobes.  Fig.  14  illustrates  this  system. 

Almost  complementary  sequences  with  positive  sidelobes 
can  be  used  to  measure  doppler.  In  Fig.  15  such  a scheme  is 

shown.  The  almost  complementary  sequences  A and  B of  length  K 

. . + 2 
give  after  summing  positive  sidelobes  at  t = - -j  K from  the 

main  lobe.  Each  of  them  is  processed  by  a matched  filter  and 

the  two  outputs  of  the  matched  filters  added.  The  output  of 

the  summer  is  connected  to  a counter  in  such  a way  that  the 

first  pulse  (main  lobe)  enables  the  counter  and  the  second 

pulse  (sidelobe)  inhibits  the  counter.  When  doppler  is 


75 


(b)  Receiving  system. 


Fig.  15.  System  to  measure  doppler. 


H 

introduced  to  the  system,  the  distance  between  the  main  lobe 
and  the  sidelobe  will  change  by  an  amount  proportional  to  the 
doppler.  The  contents  of  the  counter  is,  then,  a measure  of 
the  doppler. 

Codes  with  small  autocorrelation  sidelobes  can  be  used  in 
spread  spectrum  systems.  They  can  be  used  as  "chip"  sequences 
and  as  means  for  synchronizing  the  remote  oscillators  in  such 
systems.  They  can  also  be  used  in  surveillance  and  ranging 
systems  because  they  provide  accurate  and  unambiguous  time 
measurements . 
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VIII.  RESULTS  AND  CONCLUSION 


The  study  presented  in  this  thesis  had  as  a main  objective 
the  discovery  of  means  of  generating  "good"  codes . Good  codes 
imply  binary  sequences  having  autocorrelation  functions  with 
small  sidelobes  or  no  sidelobes  at  all. 

The  results  of  the  search  for  good  codes  in  this  study  are 
the  following: 

(1)  Almost  complementary  sequences  were  discovered. 

These  are  constructed  by  using  complementary  sequences.  The 
sum  of  the  autocorrelation  functions  has  two  sidelobes  only  of 
predictable  level,  polarity  and  position. 

(2)  Computer  programs  for  obtaining  sequences  with 
autocorrelation  functions  having  sidelobe  levels  less  than 
predetermined  desirable  levels  were  prepared  and  used.  These 
computer  programs  reveal  that  for  sequences  of  length  n = 15, 
there  are  no  codes  with  sidelobes  less  than  or  equal  to  1. 

There  is  only  one  such  code  for  n = 11  (which  is  a Barker  code) 
and  for  n = 20  three  such  codes  were  found  without  exhausting 
all  possibilities.  For  lengths  greater  than  n = 20,  an 
exhaustive  computer  search  is  impractical.  A partial  search 
can  be  made  though  by  searching  regions  of  these  sequences 
with  a special  computer  program. 

The  search  for  good  codes  resulted  also  in  the  following: 

(1)  A formula  which  gives  the  number  of  the  different 
autocorrelation  functions  for  all  possible  sequences  of  length  n. 
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(2)  A computer  program  which  can  generate  automatically 
all  sequences  of  length  n and  give  their  different  autocorrela- 
tion functions. 

Some  application  of  these  codes  were  also  considered. 

There  are  some  suggestions  for  further  research. 

(1)  Supplementary  and  cyclic  complementary  sequences 
and  their  properties  could  be  further  investigated  and  appli- 
cations developed. 

(2)  Applications  of  the  sequences  presented  here 
could  be  implemented  with  hardware. 
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APPENDIX  A 


COMPUTER  OUTPUTS  AND  PROGRAMS 


Six  computer  programs  are  included  in  this  report  and 
seven  computer  outputs  which  are  the  results  of  the  search 
for  the  "good"  codes.  In  the  programs,  the  sequence  length 
is  denoted  by  L. 

Program  1 computes  the  autocorrelation  function  of  one 
sequence,  by  comparison  of  the  elements  of  the  original 
sequence  and  its  shifted  replica.  The  length  L of  the  sequence 
and  the  sequence  (CODE  (I))  have  to  be  punched  on  separate 
cards . 

Program  2 computes  the  autocorrelation  function  of  any 
number  of  sequences  of  the  same  length.  The  number  of  the 
sequences  has  to  be  specified  in  the  loop  DO  150.  The  length 
of  the  sequences  and  each  of  the  different  sequences  have  to 
be  punched  on  separate  cards . 

Program  3 generates  all  the  possible  sequences  of  length 
L and  computes  their  autocorrelation  function.  Only  the 
length  L has  to  be  punched  on  a card.  Everything  else  is 
done  automatically.  Thus,  punching  of  the  sequences  on  cards 
is  avoided. 

Program  4 computes  only  the  different  autocorrelation 
functions  of  all  sequences  of  length  L.  This  is  done  by  the 
algorithm  explained  in  Section  VI.  The  program  produces 
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everything  automatically  and  only  the  length  L has  to  be 
specified  and  punched. 

Program  5 computes  the  different  autocorrelation  functions 
of  all  sequences  of  length  L which  have  sidelobes  less  or 
equal  to  one.  The  length  L has  to  be  specified. 

Program  6 searches  only  a region  of  all  sequences  of 
length  L and  computes  the  different  autocorrelation  functions 
with  sidelobes  less  or  equal  to  one  which  exist  in  this  region. 
Here  the  length  L and  the  starting  sequence  (SCODE  (I))  have 
to  be  specified  and  punched  on  different  cards  each.  All  the 
sequences  before  the  starting  sequence  are  ignored.  This  was 
done  because  for  long  sequences,  a complete  computer  search 
is  impractical. 

There  is  a computer  output  which  lists  the  different  auto- 
correlation functions  for  all  sequences  of  length  10 . The 
rest  of  the  outputs  list  the  "good"  codes  and  their  autocorrela- 
tion functions  for  lengths  10,  11,  12,  13,  17  and  20.  The  codes 
listed  in  all  computer  outputs  are  printed  with  the  most  signif- 
icant bit  on  the  right. 
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PRINCIPAL  VARIABLES  USED 


CODE  binary  sequence 

COR  autocorrelation  function  of  a binary  sequence 

L length  of  a binary  sequence 

ICODE  variable  used  to  complement  a binary  sequence 

JCODE  variable  used  to  reverse  a binary  sequence 

SCODE  starting  sequence  when  a region  of  the  binary 

sequence  of  length  L is  searched 
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LIST  OF  THE  DIFFERENT  AUTOCORRELATION  FUNCTIONS 


LIST  OF  ALL  DIFFERENT  AUTOCORRELATION  FUNCTIONS 
OF  ALL  SEQUENCES  OF  LENGTH  11 
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LIST  OF  THE  DIFFERENT  AUTOCORRELATION  FUNCTIONS 
OF  ALL  SECUENCES  OF  LENGTH  12 
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LIST  OF  THE  DIFFERENT  AUTOCORRELATION  FUNCTIONS 


O 

o 

o 

o 

o 

o 

o 

o 

H 

a 

«4 

o 

r-t 

1 

o 

-4 

o 

H 

o 

r-l 

o 

— t 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

z 

a 

o 

— 

a 

r4 

o 

r4  t 

o 

o 

r4 

o 

4 

o 

3 

r4 

3 

»- 

UJ 

1 

1 

m 

<_) 

< 

•4 

-j 

z 

-j 

o 

o 

o 

o 

o 

O 

o 

o 

Q 

o 

o 

o 

o 

o 

3 

< 

UJ 

UJ 

X 

3 

3 

ac 

h- 

O 

o 

CL 

o 

H 

—t 

m 

<-4 

r-t 

•4 

-4 

H 

(to 

r-t 

(to 

r-t 

(to 

r-t 

o 

UJ 

U! 

a 

1 

• 

1 

1 

z 

to 

u 

UJ 

ac 

c 

rt 

o 

o 

o 

o 

(M 

44 

o 

r-t 

o 

4* 

CM 

r-t 

o 

4t 

-i 

a 

UJ 

t- 

1 

1 

X 

3 

u. 

wo 

h 

< 

r-t 

o 

to 

ml 

«-4 

o 

r4 

O 

r-t 

o 

-4 

4 

44 

r-t 

m 

to 

1 

1 

1 

UJ 

to 

UJ 

V) 

_J 

—t 

X 

o 

o 

r-4 

o 

o 

O 

o 

o 

o 

o 

44 

M* 

o 

o 

o 

UJ 

1- 

1 

o 

to 

UJ 

z 

ut 

z 

to 

o 

-4 

•4 

H 

4) 

■-) 

r-t 

r4 

r-t 

r-t 

r-t 

*4 

cn 

-4 

UJ 

CD 

»4 

-4 

i 

3 

O 

_l 

o 

-J 

Q 

H 

o 

H 

o 

•■4 

o 

o 

CVJ 

r4 

o 

r-t 

o 

o 

(NJ 

.4 

UJ 

UJ 

0- 

Z 

1 

1 

tO 

a 

to 

o 

*— < 

a: 

u 

O 

H 

o 

H 

o 

•>4 

-4 

•-4 

o 

(to 

O 

-4 

r4 

-4 

O 

-1 

to 

-4 

UJ 

1 

— 1 

u. 

to 

<r 

X 

H 

o 

r-t 

o 

o 

OJ 

o 

CM 

-4 

o 

r-t 

3 

^t 

o 

r-t 

»- 

UJ 

UJ 

1 

1 

LL 

— » 

X 

X 

n 

»- 

►— 

O 

(to 

a 

(to 

H 

f(0 

»4 

(to 

o 

(to 

o 

fO 

o 

(to 

o 

H 

r-t 

•4 

r4 

—4 

-4 

101 


LIST  OF  ALL  (DIFFERENT  AUTOCORRELATION 
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PARTIAL  LIST  OF  THE  DIFFERENT  AUTOCORRELATION  FUNCTIONS 
OF  THE  SEQUENCES  OF  LENGTH  20 


'I 


PROGRAM  *1 


COMPUTATION  OF  THE  AUTOCORRELATION  FUNCTION 
OF  CNE  8 1 NA  T ' 


iRY  SEQUENCE 


INTEGER  COOE!  5O),COR<50> 
READ  (5  *10  )L 


RtAUin.lU  II 
10  FORMAT  (12) 


1000 


READ (5. 10 00  ) 
FORMAT! 5011  I 


(COOE ( 11,1=1,48) 


DO  40 
N=0 


1=1, L 


20 

30 

40 


100 


DO  30  J = I ,L 

IF ( CODE! J I . EQ  .CODE ! J- ( 1-1) ) ) 

N=N-1 

GO  TO  30 

N=N  +1 

CONTINUE 

COR ( I)=N 

COR ( L+l ) = 0 

MRITE!6,100>  (CODE!  I)  .1=1  ,L) 
MR  ITE(  6, 1 00  I (CORU  I.I-lfK) 
FORMAT!// ' ,30!  13  ,1X1  I 
END 


GO  TO  20 


PROGRAM  *2 

COMPUTATI ON 
OF  M BINARY 


OF  THE  AUTOCORRELATION  FUNCTION 
SEQUENCES 


(COOEI  I), 1=1,481 


INTEGER  CODE!  50),C0R(50) 

RE AO ( 5 , 10  ) L 
10  FORMAT! I 2 ) 

00  150  M= 1, 4 
READ ( 5* 1 0 00  ) 

1000  FORM AT! 50 1 1 ) 

DO  40  1*1, L 
N*0 

IF!CODEU  i^EQ  .CODE!  J-(  I-im  GO  TO  20 

K-N7i 


GO  TO  30 
20  N*N+1 
30  CONTINUE 
40  COR ( I ) *N 
CORIL+1  1*0 
K*L+i 

t*R  ITE ( 6 »l  00  ) 
WRI TE  16  ,1 00  I 
LOO  FORMAT!// 

[50  CONTINUE 
END 


(CODE!  I)  ,1  *1 » L ) 
(COR(n,I*l»K» 

• • ,30(13  , IX)  ) 


PROGRAM  #3 


AUTOMATIC  PRODUCTION  OF  ALL  SEQUENCES 
OF  LENGTH  L AND  COMPUTATION  OF 
THEIR  AUTOCORRELATION  FUNCTIONS 


INTEGER  C0DE(20>, COR (20  ) 
RE  AD ( 5 * 10  ) L 
FORMAT  ( 12  ) 

NN=2**L 

GENERATE  SEQUENCES 


00  100  II*l,NN 
DO  20  1*1,20 
20  CODE  ( I ) *0 
N I = 1 1 — 1 
DO  30  M=1,L 
JJ*L-M 

IF (NI.LT. 2**JJ) 
CODE(  JJ+1  )*1 
NI=NI-2**JJ 
30  CONTINUE 


GO  TO  30 


COMPUTE  AUTOCORRELATIONS 

DO  60  1*1, L 
N*0 

DO  50  J*I  , L 

IF ( CODE( J ) .EQ  .CODE  ( J-(  1-1 ) ) ) GO  TO  40 
N*N-i 
GO  TO  50 
40  N*N+I 
50  CONTINUE 
60  COR ( I ) = N 
COR (LH  )=  0 
K*L+1 

WR  ITE(  6,  2 00  ) (CODE  ( I),  1*1,  L) 

WRI TE(6  ,200 ) ( COR (I),I*l,K) 

LOO  CONTINUE 

>00  FORMAT(  //  ,'  ', 20(I2,2X)) 

ENO 


113 


10 

20 

30 

31 


32 


33 

34 


35 

36 


40 

50 

60 


70 

100 

200 


PROGRAM  #4 


AUTOMATIC  CALCULATION  OF  THE  OIFFER&NT 
AUTOCORRELATION  FUNCTIONS  ' 


IN  ALL  POSSIBLE 


SEQUENCES  OF  LENGTH  L 

INTEGER  C00E(20),C0R(21),IC30E(20)  ,J  CODE  (20) 
REA0(5  *10  I L 
FORMAT!  12) 

NN*2  *♦( L-l ) 


GENERATE 


SEQUENCES 
NN 


00  100  II*i,l 
03  20  I *1 » 20 
CODE  (I)»0 

ni  «ii-i 
00  30  M»1,L 
JJ*L-M 

IF  (M  .LT.2**J  J)  GO  TO  30 
CODE! JJ+1 ) *1 
NI»NI-2**JJ 
CONTINUE 

TEST  IF  SEQUENCE  ENDS  IN  0 OR  1 

IC0§i(M*={0DE(M  ) 

IF  (COOE(l).EQ.O)  GO  TO  33 

COMPLEMENT  SEQUENCE 

00  32  M*l.L 

IC  ODE(  M )*  IC  00  E (M  ) +1 

IF  (ICOOE(M)  .GT.l)  ICOOE  ( M 1=0 

CONTINUE 

REVERSE  SEQUENCE 

DO  34  M=1  »L 
JCCDElL-M+l  ) = ICOOE  ( M ) 

DO  35  M=i  « L 

IF  ( JCODE  ( L-M+ 1) -C  00  E( L-M+l ) >100,35,36 

CONTINUE 

CONTINUE 

COMPUTE  AUTOCORRELATION 
1*1  »L 


00  60 
N*0 

ificoo^uI1;  EQ.C00EU-!  I-lll  ) 

N-N-l 

GO  TO  50 

N*N*1 

CONTINUE 

COR (I) *N 

COR ( L+l ) *0 


GO  TO  AO 


K*L+i 
WRITE (6, 200) 
CONTINUE 
FORMAT  ( / ,• 
ENO 


(CODEC  I),  I *1,0, (COR (I) ,I*1,K) 
1 , 41  ( 12, IX)  ) 
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PROGRAM  *5 

COMPUTATION  OF  THE  DIFFERENT 
AUTOCORRELATION  FUNCTIONS 
OF  ALL  SEQUENCES  OF  LENGTH  L 
KITH  SIDELOBES  LESS  OR  EQUAL  TO  1 

INTEGER  CODE  <20  ),  COR  (21).  ICOOEI  20 ) , JCODE ( 20) 
REA0(5*10 )L 
10  FQRMATl 1 2 ) 

NN-2**(L-l ) 

GENERATE  SEQUENCES 

DO  100  II»i,NN 
DO  20  1-1,20 
20  CODE!  I ) =0 
NI-II-1 
00  30  M-l , L 
JJ-L-M 

IF <NI.LT.2**JJ)  GO  TO  30 
COOE(  JJ  +1  )-l 
NI-N  I-2**JJ 

30  CONTINUE 

TEST  IF  SEQUENCE  ENOS  IN  0 OR  1 

00  31  M-l  ,L 

31  ICODE(M)=CODE(MI 

IF  ( COOEil  J.EQ.OI  GO  TO  33 

COMPLEMENT  SEQUENCE 


00  32  M-l  ,L 
IC0DE(M)-IC0DE(M)+1 
IFTICODE(M)  .GT.l  ) I 
32  CONTINUE 

REVERSE  SEQUENCE 


CODE ( M )-0 


33  00  34  M*1,L 

34  JCODE(L-M+l  )-ICODE(M) 

00  35  M-l  ,L 

IF<  JCODE(L-M+1)-CODE(L-M+1)  )100,35,36 

35  CONTINUE 

36  CONTINUE 

COMPUTE  AUTOCORRELATION 

00  60  I -1 , L 
N-0 

DO  50  J-I  , L 

IF(CODE(J  J.EQ.CODEI  J-(I-l))  ) GO  TO  40 
N-N-l 
GO  TO  50 
40  N-N+l 
50  CONTINUE 

KEEP  SEQUENCES  WITH 
SI  DELOBES  LESS  OR  EQUAL  TO  l 

IF ( N .EQ.L  ) GO  TO  55 
I F ( N • GT • 1 ) GO  TO  100 
55  COR ( I ) =N 
60  COR  ( I )-N 
COR(LH)  -0 
K-L+l 

70  WRITE(6 ,200 ) (CODE!  I ) , I =1 , L ) , ( COR (I),I*1,K) 
100  CONTINUE 

200  FORMAT ( / » ' *,41(12, IX)) 

END 
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PROGRAM  *6 

COMPUTATION  OF  DIFFERENT 
AJTQCQRRE  LATI CN  FUNCTIONS  FOR  A 
REGION  OF  THE  SEQUENCES  OF  LENGTH  L 
WITH  SIDELOBES  LESS  OR  EQUAL  TO  1 

THE  SEARCH  STARTS  AT  THE 
SEQUENCE  SCODEII) , 1=1, L 

INTEGER  C0DE(20) ,COR(21I , ICOJEI 20 ) t JCODE (2 01 
READ!  5t  10  )L 

10  FORMAT! 12) 

READ  STARTING  SEOUENCE 

READ ( 5»  11 ) (SCODE (I) »I *1»  L) 

11  FORMAT  (20111 
I START*0 

DO  12  I»l,L 

12  ISTART*ISTART+SCODE(I )*2**( 1-11 
NN  »2**( L— 11—1 

GENERATE  SEQUENCES 


20  COD 
NI  i 


DO  100  II»ISTART,NN 

§?if-W° 

. II 

DO  30  M»1,L 
JJ*L-M 

IF(NI.LT.2**JJ)  GO  TO  30 
CODE!  JJ+l  1*1 
NI=NI-2**JJ 

30  CONTINUE 

TEST  IF  SEQUENCE  ENDS  I N 0 OR  1 

DO  31  M»1 » L 

31  ICODE!M)  = CODE  (Ml 

IF!  C0DE(1  » .EQ  .0)  GO  TO  33 

COMPLEMENT  SEQUENCE 

DO  32  M=I  ,L 
ICOOE!  M )a  ICODEIMM- 


IF! ICODE! M ) .GT  .1 ) 
32  CONTINUE 


1c 


ODE  (M)*0 


REVERSE  SEQUENCE 

>ICODE  ( M) 


33  DO  34  M«1,L 

34  JCOOE(L-M  + i)’ 

00  35  M*1  • L 
IF  ( JCODE!  L-M+II-COOE  (L-M+l  l 1100*35  »35 

35  CONTINUE 

36  CONTINUE 

COMPUTE  AUTOCORRELATION 
I *1  »L 
DO  50  J«I .L 

IF ( CODE ( J ».EQ.C0DE(J-(I-1)) ) GO  TO  40 
50 


00  60 
N*0 


N*  N-l 
GO  TO 
40  N*N+1 
50  CONTINUE 


KEEP  AUTOCORRELATIONS 
W IT  F SIDELOBES  LESS  OR  EQUAL  TO  1 


SCODE! 20 ) 
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K=L  + 1 

70  WR ITE ( 6j2  00 ) 
1 00  CONTINUE 
200  FORMAT  (/,'  • 
END 


(C 00 E (I)fIsL*L)  » 
41(12, IX)) 


( COR ( I ) » I sl*  K ) 
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